Isolation of microsatellite markers for the endemic Phlomis lychnitis (Lamiaceae).

BACKGROUND: Phlomis lychnitis is a mostly endemic species of the Iberian Peninsula that frequently hybridizes with the narrow endemic P. crinita subsp. malacitana in southern Spain. Whenever they coexist they form homoploid hybrid zones. Unlike hybridization at the polyploid level, the process of hybridization at the homoploid level is much less well known. METHODS AND RESULTS: In this study we report the development of 22 microsatellite markers through next-generation sequencing technologies for P. lychnitis. We characterize the genetic diversity for two populations of this species for the 10 markers that resulted to be polymorphic. Further, we check the transferability of these polymorphic markers to one population of P. crinita subsp. malacitana to verify the potential use of these markers for hybridization studies. The values of expected heterozygosity for P. lychnitis were higher than in P. crinita subsp. malacitana, and the three analyzed populations displayed negative values for the inbreeding coefficient which is compatible with the frequent instances of hybridization and introgression between species. CONCLUSIONS: This set of polymorphic markers are useful for further studies aiming at a deeper understanding of the homoploid hybrid process between these species. Additionally, this is the first panel of microsatellite markers developed for the genus Phlomis, a genus very rich in endemic species and with medicinal properties that could benefit from the use of these new markers.

De novo assembly of Phlomis purpurea after challenging with Phytophthora cinnamomi.

BACKGROUND: Phlomis plants are a source of biological active substances with potential applications in the control of phytopathogens. Phlomis purpurea (Lamiaceae) is autochthonous of southern Iberian Peninsula and Morocco and was found to be resistant to Phytophthora cinnamomi. Phlomis purpurea has revealed antagonistic effect in the rhizosphere of Quercus suber and Q. ilex against P. cinnamomi. Phlomis purpurea roots produce bioactive compounds exhibiting antitumor and anti-Phytophthora activities with potential to protect susceptible plants. Although these important capacities of P. purpurea have been demonstrated, there is no transcriptomic or genomic information available in public databases that could bring insights on the genes underlying this anti-oomycete activity. RESULTS: Using Illumina technology we obtained a de novo assembly of P. purpurea transcriptome and differential transcript abundance to identify putative defence related genes in challenged versus non-challenged plants. A total of 1,272,600,000 reads from 18 cDNA libraries were merged and assembled into 215,739 transcript contigs. BLASTX alignment to Nr NCBI database identified 124,386 unique annotated transcripts (57.7%) with significant hits. Functional annotation identified 83,550 out of 124,386 unique transcripts, which were mapped to 141 pathways. 39% of unigenes were assigned GO terms. Their functions cover biological processes, cellular component and molecular functions. Genes associated with response to stimuli, cellular and primary metabolic processes, catalytic and transporter functions were among those identified. Differential transcript abundance analysis using DESeq revealed significant differences among libraries depending on post-challenge times. Comparative cyto-histological studies of P. purpurea roots challenged with P. cinnamomi zoospores and controls revealed specific morphological features (exodermal strips and epi-cuticular layer), that may provide a constitutive efficient barrier against pathogen penetration. Genes involved in cutin biosynthesis and in exodermal Casparian strips formation were up-regulated. CONCLUSIONS: The de novo assembly of transcriptome using short reads for a non-model plant, P. purpurea, revealed many unique transcripts useful for further gene expression, biological function, genomics and functional genomics studies. The data presented suggest a combination of a constitutive resistance and an increased transcriptional response from P. purpurea when challenged with the pathogen. This knowledge opens new perspectives for the understanding of defence responses underlying pathogenic oomycete/plant interaction upon challenge with P. cinnamomi.

Chemical and genetic characterization of Phlomis species and wild hybrids in Crete.

The genus Phlomis is represented in the island of Crete (Greece, Eastern Mediterranean) by three species Phlomis cretica C. Presl., Phlomis fruticosa L., the island endemic Phlomis lanata Willd. and three hybrids Phlomis x cytherea Rech.f. (P. cretica x P. fruticosa), Phlomis x commixta Rech.f. (P. cretica x P. lanata) and Phlomis x sieberi Vierh. (P. fruticosa x P. lanata). This work describes (a) the profile of hybrids and parental species concerning their volatile compounds, (b) the suitability of ribosomal nuclear (ITS region), chloroplast (trnH-psbA), and AFLP markers to identify hybrids and (c) their competence to characterize the different chemotypes of both hybrids and their parental species. The cluster analysis and PCA constructed from chemical data (volatile oils) suggest that there are three groups of taxa. Group IA includes P. cretica and P. fruticosa, group IB includes P. x cytherea, whereas group II consists of P. x commixta, P. x sieberi and P. lanata. Volatile compounds detected only in the hybrids P. x sieberi and P. x commixta correspond to the 3% of the total compounds, value that is much higher in P. x cytherea (21%). Neighbor-joining, statistical parsimony analysis and the observations drawn from ribotypes spectrum of ITS markers divided Phlomis species in two groups, P. lanata and the complex P. cretica/P. fruticosa. In contrast to the ITS region, the plastid DNA marker follows a geographically related pattern. Neighbor-Net, PCA and Bayesian assignment analysis performed for AFLP markers separated the genotypes into three groups corresponding to populations of P. cretica, P. fruticosa, and P. lanata, respectively, while populations of P. x commixta, P. x cytherea, and P. x sieberi presented admixed ancestry. Most of the P. x cytherea samples were identified as F1 hybrids by Bayesian assignment test, while those of P. x commixta and P. x sieberi were identified as F2 hybrids. Overall, high chemical differentiation is revealed in one of the three hybrids, which is likely related with niche variation. Moreover, molecular markers show potential to identify Phlomis taxa.

Estimation of phylogenetic relationships of Phlomis species based on seed protein polymorphism

In this study, phylogenetic relationships among 39 Phlomis taxa were investigated based on seed protein profiles produced by sodium dodecyl sulphate polyacrylamide gel electrophoresis (SDS-PAGE). A total of 21 polypeptide bands were scored, of which, 19 were polymorphic among the taxa of the genus Phlomis. A distance matrix was generated from the similarity matrix which was computed by using Jaccard’s similarity coefficients, based on polymorphic bands and then an UPGMA tree was established through cluster analysis performed on the distance matrix. Genetic distances ranged from 0.00 to 0.50 within subsection Dendrophlomis; from 0.00 to 0.625 within subsection Gymnophlomis and from 0.00 to 0.769 within subsection Oxyphlomis. The UPGMA tree formed four groups. The topology of the tree is in agreement with the taxonomic view regarding the section Phlomis as it is divided into three subsections as Dendrophlomis, Gymnophlomis and Oxyphlomis based on morphological characters. The grouping pattern of the tree also indicated that subsection Dendrophlomis is more closely related to subsection Gymnophlomis than subsection Oxyphlomis.

Phylogenetic relationships among species of the subsection Dendrophlomis Bentham

This study used randomly amplified polymorphic DNA markers to determine genetic relationships among species of the subsection Dendrophlomis. Twenty accessions of the eleven Phlomis taxa were evaluated to determine genetic variability using fourteen ten mer primers selected from a 125 random oligonucleotide set. These 14 selected primers generated 85 RAPD bands that ranged in size from 200 to 1200 base pairs. Of the total bands, 88 percent (75) were polymorphic among the samples. Genetic distances among accessions were computed to produce a dendrogram based on UPGMA. Genetic distances ranged from 0.133 (between P. amanica and P. monocephala) to 0.494 (between P. chimerae and P. lunariifolia). The UPGMA tree based on distances has two major groups. The first comprised 9 taxa that were clustered into two subgroups. The first subgroup consisted of P. viscosa, P. lycia, P. amanica and P. monocephala while the second comprised P. lunariifolia, P. bourgaei, P. longifolia var. longifolia, P. grandiflora var. grandiflora and P. grandiflora var. fimbrilligera. The second group comprised 2 species, P. leucophracta and P. chimerae. Species-specific bands were observed for P. lycia, P. leucophracta, P. lunariifolia, P. bourgaei, P. chimerae and P. longifolia var longifolia.

Natural hybridization between Phlomis lycia D. Don x P. bourgaei Boiss., (Lamiaceae) revealed by RAPD markers.

Randomly Amplified Polymorphic DNA markers (RAPD) were used to assess the hybrid identity of individuals sampled as Phlomis x termessi Davis. Out of 95 primers screened, 11 primers produced reproducible amplification patterns used for discrimination of P. x termessi and their parents. Eleven primers produced 81 bands. Forty two percent of the RAPD bands existed in parents. Of the 54 bands found in P. lycia, 19 were found only in this species and 7 of these were monomorphic. Similarly, of 57 RAPD bands observed in P. bourgaei, 18 were found only in P. bourgaei and 6 of these were monomorphic. Among hybrid individuals, 35 of the 73 markers were monomorphic. Fifteen of these existed in individual parents showing that parents were homozygous for these markers. Of the 35 monomorphic bands observed among hybrid individuals, 5 were present in the samples of one of the parents and completely absent from the samples of the other; therefore, additive inheritance is indicated. Of the 5 additive bands, 1 was inherited from P. bourgaei and 4 were inherited from P. lycia. Among 38 polymorhic markers observed in hybrid individuals, 9 were new and hybrid-specific. Pollen fertility was also investigated. Mean pollen fertility for P. lycia and P. bourgaei was 93% and 97% respectively. However, mean pollen fertility for hybrids was 65% (+/-10.5).

Population genetic structure and hybridization patterns in the Mediterranean endemics Phlomis lychnitis and P. crinita (Lamiaceae).

BACKGROUND AND AIMS: The historical influence of gene flow and genetic drift after the last glacial phase of the Quaternary Period is reflected in current levels of genetic diversity and population structure of plant species. Moreover, hybridization after secondary contact might also affect population genetic diversity and structure. An assessment was made of the genetic variation and hybrid zone structure in Iberian populations of the Mediterranean Phlomis lychnitis and P. crinita, for which phylogenetic relationships are controversial, and hybridization and introgression are common. METHODS: Allozyme variation at 13 loci was analysed in 1723 individual plants sampled from 35 natural locations of P. lychnitis, P. crinita subsp. malacitana and P. crinita subsp. crinita in southern and eastern Spain. Standard genetic diversity parameters were calculated and patterns of genetic structure in each taxon were tested to fit the equilibrium between gene flow and genetic drift. Individual multilocus genotypes were subjected to Bayesian clustering analysis to estimate hybridization and introgression rates for both geographic regions. KEY RESULTS: Contrasting patterns in the distribution of genetic variation among the three taxa were found. Phlomis lychnitis showed no significant inbreeding, low genetic differentiation among populations and no evidence of isolation by distance. Phlomis crinita subsp. malacitana and P. crinita subsp. crinita showed high levels of genetic structure consistent with a pattern of gene flow-drift equilibrium. Higher instances of hybridization and introgression were detected in locations from southern Spain compared with locations from eastern Spain, matching unimodal and bimodal hybrid zones, respectively. CONCLUSIONS: High instances of historical gene flow, range expansion and altitudinal movement during the Quaternary Period, and lineage sorting can explain the diversity of patterns observed. The results suggest that P. lychnitis is the most differentiated lineage in the group; however, the relationship between the three taxa remains unclear.